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Abstract 

We study the dynamics of the reduced density matrix(RDM) of 
the field in the micromaser. The resonator is pumped by A^-atomic 
clusters of two-level atoms. At each given instant there is only one 
cluster in the cavity. We find the conditions of the independent evo- 
lution of the matrix elements of RDM belonging to a (sub)diagonal 
of the RDM, i.e. conditions of the diagonal invariance for the case 
of pumping by iV-atomic clusters. We analyze the spectrum of the 
evolution operator of the RDM and discover the existence of the qu- 
asitrapped states of the field mode. These states exist for a wide 
range of number of atoms in the cluster as well as for a broad range 
of relaxation rates. We discuss the hierarchy of dynamical processes 
in the micromaser and find an important property of the field states 
corresponding to the quasi-equilibrium: these states are close to either 
Fock states or to a superposition of the Fock states. A possibility to 
tune the distribution function of photon numbers is discussed. 
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Recent developments in the cavity electrodynamics gave rise to the cre- 
ation of a real physical device - micromaser which operates on highly excited 
Rydberg atoms pumped through a high-Q resonator 0. Existing hterature 
mostly focuses on the ideal (basic) model which is the so-called one-atom 
micromaser p, ^ ||. This device is assumed to operate in such a way 
that no more than one atom excited with the probability 1 can be found 
in the cavity at each given instant of time. The basic model is justified by 
the following assumptions: the average velocity of injection R , and the time 
of interaction r (in which a cluster passes through the resonator) are small. 
The rate of relaxation 7 of the field is low, which means a high-Q resonator. 
The coupling constant g of the field mode interacting with internal degrees 
of freedom of the atom is sufficiently large. Trajectories are assumed to be 
quasiclassical. In more exact terms these assumptions can be recapitulated 
as 



Rt < 1, r7 < 1, > 1 (1) 

The micromaser operating on a periodic sequence of iV-atomic clusters which 
are created by laser pulses in the gas of unexcited atoms is introduced in [0]. 
It is assumed that the size of the cluster is much less than the wavelength 
of the microwave radiation. Effects of finite cluster size are comparable 
in magnitude with effects of inhomogeneous field at the edges of resonator. 
The latter observation was reported in Ref. 0. In this Letter we study 
the one-cluster extension of the basic model. This means that we assume 
a point-like structure of A^-atomic cluster (i.e. the finite size effects are 
not taken into account) as well as the fulfillment of the conditions Eq.(|l]). 
This formulation generalizes greatly the basic model while leaves intact the 
simplifying assumption that process of interaction of the cluster and the field 
(within the time interval r) and the process of the field relaxation to the 
thermodynamic equilibrium (time interval T ~ 1/i?) are separated in time. 
This latter assumption allows to factorize the evolution operator of the RDM 
(see Eq. (^ below) and greatly simplifies the analysis of the properties of 
evolution operator and the dynamics of RDM. 

One-cluster model of the micromaser assumes that the A^-particle Tavis- 



Cummings Hamiltonian |Ty] 



H = Ho + V = tu (a^a + S3 + —) + 2g {^a^S. + aS+) (2) 
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is applicable. Here uj is the frequency of the quantum transition which is 
in exact resonance with the field mode. The collective spin variables S'3, 5*^ 
are the generators of the sm(2) algebra, while a\a are the creation and 
annihilation operator of the field mode, ^ = 1. 

The operation of the micromaser for each cluster is divided into two time 
intervals: the interaction time r and the relaxation time T. This means that 
the vector of the main diagonal p^^-* of RDM satisfies the following equation 



This difference equation connects the main diagonals of RDM taken at the 
instants when the /-th and (/ + l)-th clusters enter the cavity. This allows 
to understand the number of passing clusters / as a discrete "time variable". 
Here A^^e^ = -R/7, and QiNex) is the evolution operator of RDM at the relax- 
ation stage, i.e. in the empty resonator The operator iy(r) describes the 
evolution of RDM at the stage of interaction of (/+l)-th cluster with the field, 
Pat is the density matrix of A^-atomic cluster before it enters the resonator. 
The operation S'Pat means the trace with respect to atomic variables. 
In our work we consider clusters of fully excited atoms, while the field is ini- 
tially prepared in the state of the thermal equilibrium with the mean number 
of photons nf, = 0.1. In our forthcoming publication we will rigorously show 
that if there is an additive with respect to atoms and field integral of mo- 
tion [if, i^o] = and for unpolarized initial state of the cluster, then the 
dynamics of RDM is diagonally invariant. This important property of the 
evolution operator W{t) means that each (sub) diagonal of RDM in the Fock 
basis evolves independently of other elements of RD matrix. In this Letter 
we concentrate on the dynamics of the main diagonal of RDM, i.e. on the 
number of photons probability distribution function. In the space of vectors 
with components p^'-', the evolution operators Q{Nex) can be represented as 
the following matrix 



P 



(3) 





where the operator L in the matrix form reads 
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Lnm = [-{rib + i)nb{n + l)]5nm - {rib + l){n + l)Sn+l,m + nbnSn-i,m- (5) 



The matrix of the evolution operator W{t) in the Fock basis is low-triangular. 
In the present Letter we analyze this matrix by numerical methods. The 
property of diagonal invariance simplifies greatly the analysis of RDM dy- 
namics. 

The vector of the main diagonal of RDM satisfies the following difference 
equation 

P^^'^-P^'-Jr^-J^'- (6) 

Here J^'^ is the vector of the probability flux for the l-ih passage. The 
components of this vector are 

n—l n—l 

= - E E s{N)^>n"Pil + E E s{N)^,^.p'-S (7) 

n'=0 n"=n n'=n n"=0 

This vector determines the rate of change (after one passage) of the sum of 
probabilities of the photon numbers in the interval of Fock numbers between 
n = uq and n = Ui . This rate is equal to the difference of fluxes through 
the chosen boundary values, viz 

ni 

E W^^^ - p^') = 45 - (8) 

n=no 

The dependencies of the eigenvalues Wn on the number of photons are given 
in Figure 1 for the number of atoms in the cluster N=l,5,10. The interaction 
time is chosen as gr = 1.355 . The eigenvalues Wn are positive and do not 
exceed 1. Their mutual positions are defined by the parameter r and the 
number of atoms N. For the one-atom micromaser the so-called trapped 
states are known. These are the Fock states for the number of photons n 
corresponding to the eigenvalue Wn = 1 of the matrix W{t) . This number 
of photons fulfills the trapping condition 

V^^— (9) 

where x is an integer number. The trapped states do not decay in the absence 
of relaxation, and thus determine the dynamics of p'-'-' for large /. The recent 
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experimental realization of the trapped states was reported in |TT[. The 
Figure 1 shows that in the multi-atomic case there are no trapped states. 
There are however a few eigenvalues which are close to 1. The corresponding 
eigenvectors in the space of the number of photons are localized around the 
numbers n for which Wn ~ 1. Such long-living vectors is natural to call 
quasitrapped states. 

The Figure 2 shows the spectrum S{N) in ascending order. It is interesting 
to notice that the eigenvalues of the evolution operator tend to group around 
zero when the number of atoms in the cluster increases. In the hierarchy of 
dynamical processes in the micromaser the small eigenvalues are responsi- 
ble for the rapid phase of the dynamics (with respect to the discrete time / 
). The quasitrapped states corresponding to the eigenvalues in the interval 
[0.9, 1) are in turn responsible for the slow phase of dynamics. Probabilities 
of the states with corresponding photon numbers at certain stages of the field 
formation can be rather high. In Figure 2 we compare the spectrum of the 
evolution operator S{N) for the cases with {Nex = 20) and without relax- 
ation. The Figure 2 shows in particular that for bigger N the spectrum of 
S{N) is more stable towards the influence of relaxation. The relation Eq.(|^) 
describes the transition of the diagonal elements of RDM to a stationary 
state. This transition process is determined by the pumping of the cavity 
field by passing clusters as well as by the relaxation of the field. The re- 
cent literature discusses mostly 0, Q the properties of the stationary state, 
which can be achieved when a large number of clusters has gone through the 
cavity. This case corresponds to the asymptotic limit / — oo . In this work 
we concentrate on the properties of the transition process which, due to the 
existence of the quasitrapped states, are very interesting. The field rather 
rapidly "forgets" its initial state of the thermal equilibrium. The dynamics 
of the population of the Fock states shows instead the formation of long- 
living (with respect to the "time" /) quasi-equilibrium distributions. This is 
illustrated by the properties of spectrum of the evolution operator given in 
Figures 1,2. The existence of the eigenvalues close to 1 indicates consider- 
able probabilities of the Fock states with photon numbers in the vicinities of 
the maxima. The small eigenvalues correspond to the sharp depletion of the 
corresponding Fock states. The Figure 1 shows that in the chosen interval 
of Fock numbers, < n < 60, and for gr = 1.355 there are three domains 
capturing considerable probabilities. These domains, which are natural to 
call the domains of quasitrapping are localized in the vicinities of the Fock 
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numbers n = 5, 18, 40 and contain almost all the probabilities. This means 
that they are getting populated at different "moments" of "time" / in relays: 
the next domain cannot get populated until the previous one is depleted. 
This relay of populations is illustrated in Figures 3,4,5. The Figure 3 show 
for = 10 how the sums of probabilities of the Fock states change with / 
in the second (14 < n < 24) and the third (39 < n < 49) domains of qu- 
asitrapping. The rates of probability change are calculated through Eq. (|^) 
i.e. as the flux differences through the boundaries of the chosen domains. 
The Figure 3 allows to identify the following stages of the /-dynamics: a 
period of accumulation of the probability which corresponds to the positive 
values of the probability rate as well as an extended in time (/) period of 
the negative probability rates. The lasting nature of the latter period indi- 
cates that the life-time of the quasitrapped states is considerable. The rate 
of decay of the second quasitrapped state {n ^ 18) is approximately the rate 
of accumulation in the third state (n ^ 40). This means that through a 
passage of a cluster the probability is almost fully relayed from the second 
quasitrapped state to the third. Since the dynamics of the decay of the sec- 
ond quasitrapped state is slow, so is the dynamics of the accumulation in the 
third state. In Figures 4 and 5 are given dependencies on / of total popula- 
tions curves of the domains of quasitrapping. It is evident from Figures 4 and 
5 that the sum of populations of two subsequent domains of quasitrapping is 
close to 1. This again manifests the full accumulation of the probabilities in 
the domains of quasitrapping as well as the relay of probabilities indicated 
above. The Figure 6 shows the distributions of the diagonal elements of the 
RDM taken at the / -moments of maximal probabilities of the Fock states 
in the corresponding domains of quasitrapping. As follows from Figure 6 it 
is possible to govern the vector of photon number distribution. This can be 
achieved by the variation of the number of atoms passed through the res- 
onator. It is possible in particular to create states close to the Fock states 
localized at certain photon numbers. We can also report that a domain of 
the localization changes smoothly in accord with variations of parameters r, 
and N^x ■ The dependence of the stationary field on these parameters for 
A^ = 1 was discussed in Refs. 0, |^, ^ . The possibility to engineer quantum 



states is actively studied in the recent literature [|T2 
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Conclusions and discussion 



The main result of our work is the discovered possibihty to purposefully 
create in the cavity quasistable states close to Fock states. We analyzed 
the dynamics of the micromaser pumped by A^-atomic clusters 0. Our 
approach generalizes the basic model of the one atom micromaser and 
can be experimentally realized. We assumed the point-like nature of N- 
atomic clusters. This assumption can easily be realized in practice when 
clusters are created in a gas flow by focused laser pulses in the light range. 
In this case the width of the beam is of order of few microns while the 
size of the cavity can be of order of few millimeters. In our work we have 
pointed out the conditions when the time evolution of a (sub) diagonal of the 
reduced density matrix is independent of the other elements of the density 
matrix. We have investigated the properties of the spectrum of the evolution 
operator (see Figures 1,2) and discussed their connection to the properties 
of the RDM dynamics. We have discussed the hierarchy of the time scales of 
the micromaser dynamics and have shown that the sectors of the spectrum 
around zero are responsible for rapid processes while the sectors close to 1 
correspond to quasi-equilibrium. For the first time in the existing literature 
we have introduced an important notion of the quasitrapped states. The 
Figure 6 shows that these states are close to the Fock states. The domains 
in the Fock space corresponding to quasitrapping are rather narrow, their 
locations change smoothly with variations of the number of atoms in the 
cluster. This means that the overall picture of the dynamics is stable with 
respect to small variations of the number of atoms in a cluster. In our future 
work we plan to investigate this phenomenon in a greater detail as well as to 
study how the properties of the quasitrapped states depend on the choice of 
the initial density matrix of the A^-atomic cluster. 
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Figure captions 

Figure 1. The spectrum of W{0, r) for N — 



1, 5, 10 and gr — 1.355 



Figure 2. The spectrum of the evolution operator S{N) in ascending order 
for = 1, 15, Nex = 20, oo, and gr = 1.355 . 

Figure 3. The rates of change of integral probabilities of the Fock states in 
the second 14 < n < 24 and the third 39 < n < 49 quasitrapping domains 
for = 10 . 



Figure 4. Integral probabilities of the Fock states in the second 14 < n < 24 
and the third 39 < n < 49 quasitrapping domains for N = 1 . 

Figure 5. Integral probabilities of the Fock states in the second 14 < n < 24 
and the third 39 < n < 49 quasitrapping domains for = 10 . 

Figure 6. The photon number distributions of the diagonal elements of RDM 
at the / -moments of maximal probabilities of the Fock states in the second 
14 < n < 24 and the third 39 < n < 49 domains of quasitrapping. 
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